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ABSTRACT 

We present a synthesis model of the X-ray background based on the cross-correlation 
between mid-infrared and X-ray surveys, where the distribution of type 2 sources is 
assumed to follow that of luminous infrared galaxies while type 1 sources are traced by 
the observed ROSAT distribution. The best fits to both the X-ray number counts and 
background spectrum require at least some density evolution. We explore a limited 
range of parameter space for the evolutionary variables of the type 2 luminosity func- 
tion. Matching the redshift distribution to that observed in deep ^and XMM-Newton 
fields, we find weak residuals as a signature of Fe emission fr om sources in a relat ively 
peaked range of redshift. This extends the recent work of tFranceschini et all and 



emphasizes the possible correlation between obscured AGN and star-forming activity. 
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1 INTRODUCTION 

Intensive follow-up work of recent deep ^and XMM-Newton 
fields has shown that, after 40 years, the X-ray background 
(XRB) [still nolds surprises m store, wnue tfte bulk ot tne 2- 
10 keV background has been resolved into discrete sources 
(mostly active galactic nuclei; AGN), the majority of the 
power contribution to t he XRB may eme r ge at relatively 
low rcdshifts of z < 1 (Rosati et al. 2002; Hasinger 
|2002p 



z ~ 0.8 and flatten ther eafter (e.g 
[200 1|) 



Franceschini et al 



Chary fe Elbaz 



2001 



While this regime is dominated 
by obscured star-formation, reprocessed AGN accret ion ac- 
tivity may also have a non-negligible contribution. Fadda 



et al 



Brandt et al 



2002 



This is in contrast to the distribution 
of un obscured X-ray s e lected AGN which peaks at z ^ 1.5 
(e.g., Boyle et al, 1993; Vliyaji et al. 2000b), the evolution of 



QSOs observed in the optical and its dec line only at high red- 



shifts ( Boyle et al 



optic 
331 



200(1; Fan et al 



( 2002] ) found that (while the fraction of mid- infrared 
sources with X-ray counterparts is more than 10 per cent) 
at least 30 per cent of X-ray sources in the Lockman Hole 
and Hubble Deep Field-North have mid-infrared (MIR) ISO 
counterparts, and this fraction increases to 63 per cent for 
the sources detected in the 5-10 keV band. Radiative trans- 
fer modelling of the broad-band spectral energy distribution 
(SED) of a number of ISO sources has suggested that ab- 
sorbed AGN emission can indeed account for observed far- 



models which have, until now, predicted a similar peak at 
high er z (while successfully reproducing the XR B spectrum 



2001|) and XRB synthesis infra red SED peaks flWilman et al 



2000 



Franceschini et al 



2002). This evidence suggests that X-ray 'type 2' Seyferts 



Setti fe Woltjer| |l989j , |Comastri et al.| |l995j ) . Of course 



e.g. 

underlying these connections is the orientation-dependent 
unified model of AGN. 

Increasingly sophisticated models have been able to in- 
crease their goodness of fit by forcing fast evolution for ab- 



sorbed sources to lower redshift (e.g., Gilli et al 



2001 



fast evolution to 2=1.3). However, it is now likely that the 
characteristic source redshift is as low as z~0.7 (Hasinger 



2002 ; even after accounting for the effects of large scale struc- 



ture). 

On the other hand, the luminosity function (LF) of in- 
frared galaxies observed by missions such as the InfraRed 
Astronomy Satellite (IRAS) and the Infrared Space Obser- 
vatory (ISO) has been recently inferred to evolve steeply to 



and quasars may be ubiquitously detected in the 10-100 fim 
regime. 

Here we present a synthesis model for the XRB and 
X-ray number counts, based on the assumption that type 2 
obscured AGN are traced by a LF which follows the infrared 
distribution. Type 1 unobscured X-ray sources, on the other 
hand, are distributed according to the well-determined soft 
X-ray luminosity function (XLF). We include a distribution 
of source spectra, accounting for full Thomson scattering 
cross-sections and explore a range of parameter space for 
volume emissivity evolution with redshift. We comment on 
the need for various classes of sources in matching the obser- 
vational constraints, and investigate the effects of including 
an Fe Ko emission lin e in the template spectra. T his follows 
on from the work of Franceschini et al. (2002, hereafter 
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F02), who adopted a fixed fraction of a composite MIR lu- 
minosity function for the density and evolution of a single 
type 2 X-ray SED, and were successful at reproducing the 
observed peak of the redshift distribution in deep fields. Fi- 
nally, we construct far-infrared SEDs under the assumption 
of complete reprocessing of the AGN blue bump radiation to 
predict obscured AGN source counts at 70 /um, and discuss 
implications for unification models. 

For the purposes of this work, we ignore distinctions in 
the optical and X-ray classifications of AGN: i.e., type 1 = 
unobscured (or mildly obscured, see next section) and type 
2 = obscured. We assume Ho=50 km s _1 Mpc -1 , f2M=l and 
Q A =0- 



2 THE MODEL 
2.1 Spectra 

The AGN X-ray spectra are assumed to be power-laws with 
a slope fixed at a=0.9 (Fe oc E~ a ) and an exponential cut- 
off at Ec = 360 keV. A reflection hump from cold material 
at sol ar abundance is added by em ploying the PEXRAV 
code (Magdziarz & Zdziarski 1995), assuming uniform 27r 
reflection and inclination angle 9 = 60°. Iron emission lin es 
are not included here, and are discussed separately (§ |3,3| ), 
Obscuration by material of solar abundance local to the 
source is modell ed (with log7VH=20..,25; step .5) using the 
cross-sections of 



Morrison fe McCammon ( 1983 ) . The impor- 



tance of Compton down-scattering for Thoms on-depths (r 
>1 (lo gA r H>24) has bee n discussed previou sly [ belotti et al 



(1995, hereafter C95); Wilman fe Fabian ( 



199S, hereafter 



WF99)] and we model the full Klein-Nishina cross-section 
using the same monte carlo code as the latter authors, again 
assuming solar abundance. 

Finally, 2 per cent of the intrinsic continuum is added 
to all spectra to represent the component scattered into the 
line-of-sight. See Fig 1 of WF99 for an illustration of the 
constructed spectra. 

The ratio of type 2 : type 1 sources in the local Uni- 
verse is thought to be around 4 and could be larger than 10 



(Maiolino & Rieke 



1995 



Matt et al 



200C). We parametrize 



The LF of type 2 X-ray sources is completely unknown, 
due to the difficulty of compiling an unbiased and complete 
sample. As mentioned, we assume that these are traced by 
a MIR local luminosity function (LLF). A number of such 
LLFs have been d etermined , and here we adopt the one 

(1998), based on incorporation 



Ku et al 



found at 15 fim by 
of IS O-determined SEDs in to 14 06 IRAS galaxy o b serva - 

~ (pool) 



Chary & Elbaz (2001) and 



Franceschini et al 



tions. 

have used this LLF (partially or fully, in combination with 
other LLFs) to successfully reproduce the slope of the ob- 
served MIR number counts at faint fluxes with appropriate 
evolution. We have not explored the dependence on the var- 
ious LLFs, since our purpose here is only to demonstrate the 
viability of the X-ray:IR connection. 

There are then two major unknowns in this work, for 
which there are no firm observational constraints: 1) the rel- 
ative X-ray :IR intrinsic luminosity ratio for each SED; 2) the 
fraction of IR sources which are X-ray detected AGN, as a 
function of logA^H- The observed X-ray:IR ratio will depend 
on the amount of obscuration and larger absorbing columns 
of gas and dust will deplete more X-rays and reprocess larger 
amo unts of optical/UV radiation to the IR. Using the stan 
dard 



Elvis et al 



(1994) median SED for radio-quiet quasars 
(with the normalization of the LF discussed below) leads 
to an underestimation of the XRB spectrum by as much 
as a factor of 3 (or overestimation of the counts, if we cor- 
respondingly correct the LF upwards). However, there are 
suggestions that the X-ray:IR ratio can be similarly large in 
the case of Seyferts. This can also be interpreted as assum- 
ing that the bolometric correction factor (i.e., L2-10 / Lboi) is 
closer to 10 per cent, as compared to 3 per cent for quasars, 
where L2-10 is the intrinsic X-ray luminosity and Lboi is the 
bolometric luminosity. For instan ce, the general t rend of lu- 



minosity ratios for the sample of Banders et al. (1989 



this distribution as (logNu^ ■ Type 1 sources are taken to be band by Piccinotti et al 



their Fig 9) suggests this possibility. Indeed, if the XRB is 
dominated by Seyferts as opposed to quasars, then a larger 
bolometric correction may be appropriate. With this approx- 
imation, we find that L2-10 ~ iis^m- 

To determine the fraction of IR sources which will be 
X-ray detected AGN (the rest are probably obscured star- 
bursts, as discussed in the introduction), we choose to match 
the MIR LF to the X-ray LF det ermined in the 2-10 keV 



those with logA r H^22, with the rest being type 2. We note 
that this renders our model effectively 'unified', since such 
a parametrization can be expected from the angular distri- 
bution of torus column density, where the maximum optical 
depth is along the equatorial plane with the most probable 
orientation (discussed in C95). However, this qualitatively 
matches the increasing fraction of highly obscured so urces 
(at least to logA r H=24) observed by Risaliti et al. (1999) and 



also allows us to explore a range of possible distributions. 



2.2 Luminosity functions and their evolution 

The type 1 XLF that we use is the observed-frame 0.5-2 



keV double power-law form found by Miyaji et al. (2000b) 



with a density evolution which is luminosity-dependent. The 
parameters of the LF (both local and evolutio nary) that we 



adopt a re those of the BAS1 model quoted in Miyaji et al 



(2000a). The only change is the relative weighting of type 1 
sub-classes (with different Nh) as described above. 



(1982), and adopt this as the 
XLF for each Nh sub-class. We do this by simply scaling 
so that the two match at 10 44 erg s _1 . No SED correc- 
tions are needed, since the 2-10 keV and 15/xm luminosities 
are approximately equal. The adopted LLFs are shown in 
Fig While this prescription may be rather crude, it is self- 
consistent in the sense that the local X-ray space density 
(from all Nh sub-classes) adds to a few per c ent of the IR 
space-density, as found by Fadda et al. (2002, discussed in 
the following sections). 

For the evolution with redshift of the type 2 LLF, we 
investigate power-law emissivity evolution [as (l + z) p ] out to 
Zcut and a constant weighting thereafter to z mBX (set equal 
to the weight at z cu t)- 



2.3 Model fitting 

Integration has been carrie d out from logL m jn =41.7 to 

L refers 



logL max =47.5 (the limits of Miyaji et al 



2000b 



to the observed 0.5-2 keV band for type 1 sources and the 
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2-iO keV (Piccinotti et al. 82) 
■ 15 iim LLF rcsealed at 10 44 erg s 1 
Miyaji et al. (2000) 0.3-3.5 keV LLF 



org s 1 Mp c , in agreement with the values derived by 



0.100 1.000 
L. M (10" erg S -') 



Figure 1. The local luminosity functi ons relevan t to th is paper. 
The dotted line is the 15 /im LLF of Xu ct al.| <199$i. This is 
normalised to the 2-10 keV X-ray LLF of |pkcmotti ct aL ( 1982[ 



dot-dashed line) at 10 erg s — the result is the solid line, 
which is the adopted type 2 LLF for each Nn sub-class without 
the appropriate weighting factors. The combined LLF for all the 
type 2 sources will depend on the type 2 : type 1 ratio, which will 
vary with (3. The dashed line is the ad opted LLF of t ype 1 s ources 
in the observed 0.5—2 keV band from Miyaji et al.| ( l2000b| ). 



rest-frame 2-10 keV band for type 2 sources in units of erg 
s _1 ) for all the source sub-classes, with the variables being 
the type 2 parameters: z max , z cut and p. In accordance with 
C95 and WF99, we normalize the generated spectrum by 
a factor (determined during fitting) between 0.7 and 1.5, in 
order to account for uncertainties in the absolute level of the 
background, deviations from the 'true' luminosity function 
and degeneracy of cosmology. A good fit is assumed if the 
deviation of the normalized spectrum on a grid of energies 
at 10, 20, 40, 50, 75 and 100 keV is ^10 per cent, while a 
maximum contribution of 97 per cent at 5 keV is required, 
in order to allow for contributions due to clusters, which we 
do not explicitly model here. 

We further discuss a model that provides a good de- 
scription of the X-ray number counts and redshift distribu- 
tion. However, an explicit fit is not performed due to the 
number of parameters involved. 



Miyaji et al. (1994) by cross-correlating the XRB with IRAS 



galaxies. 

More generally, the region of parameter space explored 
is shown in Fig H. For the models shown, z ma x values of 
0.5... 5 were explored. The shaded regions show models which 
deviate fr om the analyt i c fit t o the observed XRB deter- 
mined by Glruber et al. (1999) by at most 12 per cent at 



the energy values mentioned. The smaller black squares lo- 
cate models where the maximum deviation is 10 per cent. 
The range of acceptable models is relatively narrow, since 
advancing along either axis corresponds to increasing source 
weighting. WF99 pointed out that parameter space could be 
opened up by shifting the individual spectral peaks above 30 
keV by allowing super-solar metallicities. Including the ef- 
fects of coronal Comptoniza tion of the reflection hump may 
also help to a small extent (Petrucci et al. 2001). 

The other observational constraint that we match is 
the integral number counts in the 2-10 keV band. We have 
explored splitting emissivity evolution into luminosity and 
density evolution (with power-law indices pl and p„ respec- 
tively) without affecting the XRB spectrum itself. While 
there is degeneracy between these possibilities, we find that 
some density evolution is required to delay the flattening of 
the integral number counts to below 10 -16 erg s _1 cm -2 
(as also pointed out by F02). Models with pure luminosity 
evolution produce a much flatter slope in the counts at faint 
fluxes. We present one model with pl~2 and p„=2 in Fig ^. 



The observed counts are from Rosati et al, (2002) 



The redshift distribution of sources predicted by this 
model is shown in Fig ^, where the histogram is the (un- 
normalised) source distribution from Eiasinger (2002, their 



Fig 6). The spectroscopic completeness fraction reported by 
Hasinger is about 60 per cent - the rest are likely to be type 
2 sources at intermediate to high redshifts with weak (or ab- 
sent) emission lines. Low redshift sources are likely to have 
a high identification completeness. We have thus normal- 
ized the observed distribution to the model distribution at 
z = 0.5 (avoiding z — 0.7 due to the presence of large scale 
structure). The fraction of identified sources beyond z = 1 is 
about 40 per cent of the total predicted number. Exact com- 
parison is difficult due to the mixture of flux limits of the 
surveys from which the identified sources are drawn. As also 
mentioned by F02, at least part of the shortfall of sources 
at z < 0.5 could be due to the contribution by clusters. 



3 RESULTS 

3.1 Matching observational constraints 

The fit to the XRB spectrum for one of our 'good' models 
is shown in Fig |^, for which z cut =0.7, z max =4.5 and p—i, 
with a normalization factor of 1.44. The contribution from 
various classes of AGN is shown as well, with type 1 sources 
dominating below 2 keV, while Compton-thin type 2 sources 
dominate beyond. Compton-thick sources with log7VH>24 
never dominate the XRB due to the effect of Thomson down- 
scattering. 

The local number density as a function of logiVH re- 
quired by this model is shown in Fig |^. The weighting pa- 
rameter (/3) = 5, implying a local type 2 : type 1 ratio of 5.6. 
The local 2-10 keV volume emissivity implied is 4.5 xlO 38 



3.2 The need for highly obscured sources 

The major discrepancy between the observed and the mod- 
elled XRB spectrum occurs between 10 and 20 keV, close 
to the position of the redshifted peak (from 30 keV) of a 
source at the characteristic redshift of 0.7 (Fig|2|). The max- 
imum deviation at 14 keV is slightly more than 10 per cent. 
Although there is no doubt of the need for obscured ac- 
cretion in producing the hard spectrum of the XRB, the 
relatively sharp peak of the Compton-thick contribution 
(logA r H>24) at ~20 keV (dotted line) suggests that remov- 
ing these sources may decrease the maximum percentage de- 
viation of the model. We have tried setting logiV H nax =24 and, 
indeed, find many solutions which fit the XRB spectrum 
well. However, removal of these sources effectively makes 
the total spectrum softer, and thus forces a much lower 
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Figure 2. The crosses represent the observed XRB spectrum from Grubcr et al. (1999). Our model has been fit to the analytic form from 
the same paper between ~10 and 100 keV. The solid line is the total spectrum with a normalization factor of 1.44. Contributions from 
type 1 sources (dashed), type 2 Compton-thin sources with 22:<Clog7VH<24 (dot-dot-dot-dashed) and Compton-thick sources (dotted) are 
shown. 



Zmax (with appropriate renormalization) to reproduce the 
flat slope below ~30 keV. This leads to number counts be- 
ing severely underestimated. In addition, the contribution of 
Compton-thick sources is important to the 2-10 keV source 
counts at the very faintest fluxes (Fig ^), and this should 
be even more so for counts at harder energies. Thus, they 
cannot be easily excluded from the model. 

We also note that marginally acceptable models with no 
type 2 quasars (i.e. sources with logA r H>22 and logL>44) 
can be found. Their contribution is relatively small (due to 
the decrease in their number density by >10 4 compared to 
logL=42; Fig|l|), mainly showing up in the number counts 
at brighter fluxes above 10 -15 erg s _1 cm -2 . 



3.3 Features in the spectrum 

One consequence of having a 'characteristic' source redshift 
around which the bulk of the background is produced is that 
an imprint may be left on the background spectrum itself. 
No such characteristic features (incl uding any Fe emission 



lines) have been found to date (e.g., Schwartz 



1992 



Fabian 



& Barcons 1992 and references therein) and this has been 



taken as strong evidence for a truly cosmological popula- 
tion of sources which is well spread-out in redshift, thereby 
smearing out any features. We investigate the effect of incor- 
porating an Fe Kq emission line with our relatively-peaked 
redshift distribution. 

The lines are assumed to be narrow gaussian profiles 
with a fixed width a= 10 eV (In practice, the resolution of 
the energy grid was such that the lines had a relatively sharp 
peak). The rest-frame equivalent-widths (EW) of lines in 
type 1 sources has been taken to be 150 eV ( George fc Fabian 



L_ 
no 



Figure 3. Local space density of sources implied by the model in 
Fig |^. The light, thin line shows the differential density per unit 
\ogNn. The solid line is the cumulat ive density. For c ompa rison, 
the number density predicted by the Piccinotti et al. (1982) XLF 
(integrated over the luminosity limits stated therein) is shown as 
the dashed line. 



1991). For type 2 sourc es with logA f H>23.25, we match 
to the relation found by Leahy & Creighton ( 1993[ ) where 



the EW is 180 and 900 eV for logiV H =23.75 and 24.25 re- 
spectively, while for sources with logA r H=22.25...23.25, they 
found an EW less than 150 eV. Thus, for these sources, we 
fix EW=150 eV as well. For the highest column densities, 
Leahy & Creighton found column-densities greater than 10 
keV. To be conservative, we choose to fix the Fe Ka EW to 



X-ray background synthesis: the infrared connection 5 



10000 E 



1000 



100 - 



A 




10 L ' J 10 10 

S (2 - 10 keV) [erg s"'cm" ! ] 



12 



Figure 5. 2-10 keV source counts predicted by the model of Fig ^ (s olid line), wit h cont ributions from various source types similarly 
shown. The solid circles are a sample of the source counts reported by Rosati ct al,| (J200 2j). Three repr esentative error b ars (with hats) 
are s hown. The diamonds with error-bars are from (with increasing flux) ASCA dOgasaka ct"aT] [iggg ) and BeppoSAX (]g iommi et al 
2000| ). 



1000 



100 



10 




Figure 6. Source distribution as a function of redshift, calculated in steps of 0.2 in z (solid line) for a flux limit S(2-10 keV)>5xlO 16 
erg s _1 cm -2 . The contribution o f type 1 s o urces is shown as the dotted line, while that of all type 2 sources is shown as the dot-dashed 



line. The observed distribution of 
(see text). 



Hasingei (2002, histogram and black circles) has been normalised to the predicted number at z = 0.5 
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2 4 6 8 10 



f 

Figure 4. The range of acceptable models which fit the XRB 
spectrum only. The shades of grey represent different values of JVh 
weighting index /3. Lightest: /3=2; Intermediate: /3=5; Darkest: 
[3=8. These represent deviations from the observed spectrum of 
less than 12 per cent at 10, 20, 40, 50, 75 and 100 keV, after 
normalization as described in the text. The filled black squares 
represent models with deviations of less than 10 per cent. z max 
values for these models typically range from 3—5. 
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Figure 7. Percentage difference in the modelled spectrum when 
Fe Kct lines are included with conservative assumptions (see text), 
compared to the best model without the emission lines. The peak 
(and bulk of this feature) is due to type 2 sources at a character- 
istic redshift of 0.7, as described. The much weaker second peak 
corresponds to the break at 2=1.5 of type 1 sources. The sharp- 
ness of the peaks is partly due to the fact that a low-resolution 
grid was used for modelling the lines. 



6 keV for logA r H=24.75 (this was the maximum EW found 



Levenson et al 



2005). This 



in a sample of Seyfert 2s by 
is the strength of the line generated from fluorescence in 
the obscuring matter itself. Fig ^ shows the effect produced. 
Compared to the best model without the Fe line included, 
we find a maximum difference of more than 3 per cent close 
to 6.4/(1 + z)=3.8 keV, where 2=0.7 is the characteristic 
source redshift. Though inclusion of cluster emission may 
dilute the EW of this feature on the XRB spectrum itself, 



especia 



, ly at soflei eneigies, we emphasize thai such a weak 
feature ' should be hed for, esp e cially if th e (observed) 

redshift distribution peaking at 2=0.7 is to be believed. 



3.4 Obscured AGN contribution to IR counts 

Finally, we make use of the correlation between X-rays and 
the IR to make predictions of type 2 AGN source counts at 
70 fim (the MIPS instrument aboard the forthcoming SIRTF 
mission has a 70 [im 5a sensi tivity of 1.3 mjy in 500 s; e.g., 
Brandl fe the IRS team] |200C| ) . 

For fc-corrections, we constr uct models with t he public 
radiative transfer co de DUSTY ( Ivezic et al. 199!:), with all 



the assumptions of Wilman et al. (2000). The main vari- 
able values assumed are : r out /ri n =500, the relative scales 
of the spherical dust distribution; Ti n =1000 K, the dust 
temperature at r ln and t(0.3 fim)=10... 10000, the optical- 
depth. These then correspond to type 2 X-ray sources with 
logA r H~22...25, assuming a Galactic dust:gas ratio and so- 
lar abundances. DUSTY cannot currently compute a ra- 
diative transfer solution through a torroidal dust geome- 
try. Although we are thus biased towards fainter sources 
(and our estimates will be lower-limits), we note that new 
observations are finding many objects which have a higher 
covering fraction than that inferred from a torus. Absence 
(or weakness) of optical and near-infrared emission lines 



in many newly-discovered ^and XMM-Newton sources sug- 
gests a dust covering fraction large enough to completely 
obscure each source. Indeed, the flat shape of the observed 
XRB spectrum suggests a very large covering fraction (e.g., 



Fabian fe Iwasawaj 199E deduce this to be approximately 



85 per cent). Moreover, good fits to the broad-band SEDs of 
optically-red and X-ray hard sources have been obtained un- 
der the assumption of c omplete obscuration a nd / or very 



low scatter ing fraction (Wilman et al 



et al 



200C 



Franceschini 



2002). The case of spherical obscuration may be valid 



for a non-negligible fraction of distant type 2 AGN. We re- 
turn to this point in the next section. 

To determine the IR:X-ray ratio, we determine the ratio 
of the bolometric to 70 /im (monochromatic) luminosity for 
each SED, and convert this t o an intrinsi c 2-10 keV luminos- 
ity using the ratio found by Elvis et al. (1994) for a sample 



of local, unobscured, radio-quiet quasars. We can now use 
the type 2 X-ray LLF and the evolution derived for our il- 
lustrative model and predict 70 /im source counts above a 
limiting flux of 0.1 mjy. These are sho wn in Fig | g| Com- 
paring this with the prediction made by Franceschini et al 



( pOOl , their Fig 20), we find that at the faintest fluxes, type 
2 AGN will constitute less than 1 per cent of the cumulative 
IR counts. However, our estimate must be treated as a lower 
limit since we have not modelled starburst contributions to 
the SED (which would make the fc-correction higher - es- 
pecially the most obscured sources - and increase number 
counts significantly). 

If, on the other hand, the dust:gas ratio is lower than 
the Galactic va l ue (a s has been seen at low redshift; e.g.. 



Maiolino et al 



2001), the IR counts will increase signifi- 



cantly. For instance, assuming for simplicity that the typi- 
cal extin ction due to dust is uniformly lower by a factor of 
10 (e.g., Granato et al. 1997) than that expected from the 



observed amount of X-ray obscuration (and if the dust:gas 
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Figure 8. (Top) DUSTY spectra with the parameters described 
in the text and X-ray logiVH=22 (dotted), 23 (dashed), 24 (dot- 
dashed), 25 (dot-dot-dot-dashed), corresponding approximately 
to no. 3;zm=10... 10000. (Bottom) Type 2 AGN source contribution 
to the 70 fim source counts with the fc-corrections of the spectra 
in t he t op figure, and density distribution of the model described 
in 5 3.1 Total counts - solid line; individual contributions - as in 



the top figure. 



ratio is Galactic), we obtain number counts which are higher 



by a factor ol two at the lamtest flux level. 



4 DISCUSSION 

We have presented a model in which a distinction is made 
between type 1 and type 2 X-ray source distribution and 
evolution, and type 2 sources follow a (rescaled) luminosity 
function of MIR sources. With simple power-law emissiv- 
ity evolution, we have obtained relatively good fits to the 
XRB spectrum over 5~100 keV. Some density evolution is 
required to match the increasing integral source counts at 
faint fluxes. We have not fit parameters relevant to type 
1 sources, although their relative fraction does vary with 
changing (3. 3-5 per cent of IR sources are found to be ob- 
scured AGN locally (4 per cent for the model in Fig ^ and 
this increases to 7-13 per cent (of the local IR number den- 
sity) at z = 0.7, depending on the value of j3. We note that 



the illustrated model is not the 'best' one in the sense of 
having the minimum deviations from the XRB spectrum 
of all possible models. However, this does have moderate 
success in accounting for all observational constraints (the 
XRB spectrum, X-ray number counts and observed peak 
of the redshift distribution). The main uncertainties remain 
the exact LF of type 2 sources and the IR:X-ray luminosity 
ratio for each Ah sub-class. 

We find residuals on the XRB spectrum due to intrinsic 
Fe emission which are not significantly smeared out by the 
peaked redshift distribution. This feature is weak and at the 
level of a few per cent. It lies close to t he low er range (3 keV) 
of energies over which Marshall et al. (198C ) found no signifi- 



cant residuals; its non-detection, therefore, is not necessarily 
in conflict with our model, and further effort should be made 
to either confirm or rule out this feature observationally. 

We have extended the work of F02 by using a distribu- 
tion of absorbed X-ray spectra with full Thomson scattering 
cross-sections. In so doing, however, we have less certainty 
on the IR:X-ray luminosity ratio for each Ah sub-class. This 
could explain our lower estimate (4 per cent) of the local 
fraction of IR sources which are obscured AGN. We have 
investigated a range of parameter space for the type 2 evolu- 
tionary parameters, and commented on the need for various 
classes of sources. Differences with F02 in the determined 
redshift distributions are partly due the different IR LFs 
used as a starting point and to differences in cosmology. For 
instance, the number of type 2 sources per deg 2 in a non- 
zero lambda cosmology with Qm = 0.3 and Qa = 0.7 is lower 
by about 15 per cent at z = 4 (refer to Fig ^j; this effect is 
larger at high redshifts due to the rapidly decreasing comov- 
ing volume). Additionally, adopting a higher flux threshold 
of 8 x 10~ 16 [cgs] leads to 40 per cent less counts at z = 4 
(the limiting fluxes of the observed redshift distribution vary 
between the samples, which makes exact comparison diffi- 
cult). 

The intimate connection between star-formation and 
AGN accretion activity has been known for some time now, 
especially in t he m ost massive galactic system s (see e.g., 



Genzel et al 



1998 



and 



Magorrian et al 



1996 for differ- 



ent lines of evidence). Phenomenological models which self- 
consistently account for the symbiotic gro wth of a black hole 
and its host bulge have been explored (e.g., ^abian 



Rees 



lOOC; Silk 



1998) and can reproduce observed correlations such 



as t he mass-dispersion r elation of Gebhardt e t al. (200C) 

( pOOq ) have 



Merritt & Ferrarese (2001b) 



Haiman et al 



and 

shown that, in a primordial cloud, enhanced production of 
photoelectrons due to emission of hard X-rays from a mini- 
quasar facilitates the formation of H~ , leading to increased 
abundance of H2 which allows cooling and collapse of dense 
virialized clouds, possibly leading to star-formation. This ef- 
fect may be stronger for halos with higher metallicity. 

What the new deep field observations suggest, is a uni- 
versal peak phase of accretion activity at z ~ 0.7, which is 
delayed from the phase of unobscured quasar activity peak- 
ing at or beyond z ~ 1 — 2. Recent determination of the star- 
f ormation history also shows a flattening beyond z ~ 0.7 
( |Chary fe Elbaz ), and it has already been suggested that 



an increased rate of mergers at these intermediate redshifts 
could supply additional gas to trigger as well as feed nuclear 
star-formation and also lead to obscured accretion. While 
an increase in the merger rate has been observed to 2=1 (|Lej 



8 P. Gandhi and A.C. Fabian 



Fevre et al. 2000), a detailed study of this merger rate would 



be essential to accurately determine the higher redshift c on- 
nectio n with feeding obscured AGN. We also note that Xu 
(2000) found a higher z cut = 1.5, based on the direct (rather 
than derived) 15 /im ISO LF. Follow-up observations of dis- 
tant SIRTF sources may resolve this issue within the next 
few years. 

The local mass density in black holes based on the 
redshift distribution of our model is found to be pbh = 
7.6xl0 5 (j7/0.1) _1 (//0.03) _1 M Mpc" 3 , predicting a bolo- 
metric accretion energy density of 1.3 x 10 58 (//0.03) _1 erg 
Mpc -3 (77 is the accretion efficiency and / is the 2-10 keV 
: bolometric correctio n factor adopted from the radio-quiet 



Elvis et al. 



quasar compilation of 

This is more than a factor of 2 larger than the mass 



1994| ) 



density found by Fabian & Iwasawa (199£) from simple ab- 



sorption correction of the observed XRB spectrum (adopting 
a median redshift of 0.7), and less by the same fac tor as com- 



pared to the upper range of densities quoted by Elvis et al 
(2002), based on better estimates of the bolometric correc- 
tion, < z >=2 and assuming 77=0.1 (We compare with the 
upper limit since our estimate would al so in crease with the 
new c orrectio ns) . Yu & Tremainc (2002) and Merritt & Fer- 
rarese ( 2001a ) find densities which are lower by a factor of 2- 



4 (after homogenizing the adopted cosmology) base d on the 



Yu & Tremainc 



conclude 



locally observed Mbh — relation 
that efficient accretion in optically-bright QSOs can account 
for the calculated pbh- In this regard, we note that a signifi- 
cant fraction of our mass density is due to obscured accretion 
at high redshifts (we integrate up to z m ax=4.5), and if we 
adopt a Zmax=l-5 for type 2 sources, pbh drops to 4.9 x 10 5 
Mq Mpc -3 without affecting the XRB spectrum by the same 
factor. However, the flat slope of the XRB spectrum unam- 
biguously suggests the presence of powerful objects which 
are highly obscured, and this is difficult to reconcile with 
the spectra of type 1 QSOs. Again, the results from SIRTF 
surve ys (e.g ., Btorric-Lombardi fc The FIRST Look Survey 



Team flieke & The Sirtf Guaranteed Time Observers 



2001) should prov i de im portant constraints. 

Fabian et al. ( 1998 ) discussed the possibility that nu- 
clear star-formation could provide the energy required to 
keep a circum-nuclear torus puffed up with the large covering 
fraction implied by the XRB spectrum. However, the com- 
plete absence of emission lines (from either star-formation 
or AGN activity) in b oth the optical (e.g. , Hasingei 2002) 



and the near-infrared (Gandhi et al. 2002) in many X-ray 
obscured sources suggests that star-formation is either ab- 
sent, or highly reddened. Although larger samples of ^and 
XMM-Newton follow-up spectra are required for quantita- 
tive analysis, it seems clear that obscuration at high redshift 
must be orientation-independent and that the simple AGN 
unification based on orientation cannot hold at increasing 
redshift. Whether the difference between distant type 1 and 
type 2 AGN is due to nature (the conditions during forma- 
tion) or nurture (effects of the environment, e.g., mergers) 
or a mixture of these effects, needs to be understood. 
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